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In applying vibration-transit (V-T) theory of liquid dynamics to the thermodynamic properties 
of monatomic liquids, the point has been reached where an improved model is needed for the small 
(~ 10%) transit contribution. Toward this goal, an analysis of the available high-temperature exper- 
imental entropy data for elemental liquids was recently completed [D. C. Wallace, E. D. Chisolm, 
and N. Bock, Phys. Rev. E 79, 051201 (2009)]. This analysis yields a common curve of transit 
entropy vs. T/Otr, where T is temperature and 9 tr is a scaling temperature for each element. In 
the present paper, a statistical mechanics model is constructed for the transit partition function, 
and is calibrated to the experimental transit entropy curve. The model has two scalar parameters, 
and captures the temperature scaling of experiment. The calibrated model fits the experimental 
liquid entropy to high accuracy at all temperatures. With no additional parameters, the model 
also agrees with both experiment and molecular dynamics for the internal energy vs. T for Na. 
With the calibrated transit model, V-T theory provides equations subject to ab initio evaluation for 
thermodynamic properties of monatomic liquids. This will allow the range of applicability of the 
theory, and its overall accuracy, to be determined. More generally, the hypothesis of V-T theory, 
which divides the many-atom potential energy valleys into random and symmetric classes, can also 
be tested for its application beyond monatomic systems. 

PACS numbers: 65.20.De, 05.20.Jj, 05.70.Ce, 61.20.Ne 



I. INTRODUCTION 

Our physical understanding of the motion of atoms in 
real liquids is notably deficient compared to lattice dy- 
namics theory for real crystals. The present work stands 
in a long line of research aimed at improving the liquid 
theory. We shall consider only monatomic liquids, and 
equilibrium thermodynamic properties, though progress 
has recently been made in nonequilibrium theory as well 

To rationalize experimental data, theory must be based 
on a physically realistic interatomic potential. Pseudopo- 
tential perturbation theory provides such interatomic po- 
tentials for the nearly- free-electron metals Q . These po- 
tentials have been valuable in developing liquid dynamics 
theory, because the metals to which they apply are very 
well studied experimentally. In the early days, pseudopo- 
tential perturbation theory was developed by comparing 
theory with experimental data for crystals. Phonon dis- 
persion curves were calculated for Na by Sham (ij, and 
for Al by Harrison Q . Pseudopotential form factors were 
compared with Fermi surface data by Ashcroft [5|], and 
compressibilities and binding energies were calculated by 
Ashcroft and Langreth [6J. Then, using a well-tested 
pseudopotential for Rb, Rahman demonstrated the abil- 
ity of molecular dynamics (MD) to produce accurate re- 
sults for statistical mechanical properties of the liquid 
0, A host of MD calculations followed, showing 
good agreement with experimental data for nearly-free- 
electron liquid metals. These included the structure fac- 
tor of alkali metals [1, [l(| , thermodynamic properties of 
Na [111, a nd structural and thermodynamic properties of 
alloys [|3. 

The next step in liquid dynamics theory for real ma- 



terials was the introduction of ab initio MD by Car and 
Parrinello [la ]. This is still based on pseudopotentials 
in electronic structure theory, but the pseudopotential is 
treated numerically instead of as a perturbation. The 
early development of computational methods is reviewed 
by Payne et al. [l4| • Ultimately ab initio MD calculations 
of the nuclei moving on the ground state adiabatic po- 
tential surface have produced results in good agreement 
with experiment for the melting properties of Al (de Wijs 
et al. 15]), structural and dynamic properties of liquid 
Fe under earth's core conditions (Alfe et al. [l6|), and 
the pair distribution functions for groups IIIB-VIB ele- 
mental liquids (Kresse [l7j : see also Chai et al. [l8| for 
Ge). An extensive study of Na at high compression re- 
veals significant electronic structure changes in the crys- 
tal (Neaton and Ashcroft [l!|) and in the liquid (Raty et 
al. [20(), and a change from normal to anomalous melting 



(Gregoryanz et al. 21]). Finally, ab initio Monte Carlo 



calculation of liquid free energy is also being developed 
(Greeff and Lizarraga [H , Greeff H|). 

The goal of V-T theory is to develop a Hamiltonian 
formulation capable of analyzing the motion of atoms 
in the liquid state. The theory is based on a funda- 
mental hypothesis, and an advantageous decomposition 
of the atomic motion. These propositions are formu- 
lated to rationalize a larg e body of experimental data 
for monatomic liquids [24. [2o| . The hypothesis, a sym- 
metry classification of potential energy valleys, was ar- 
gued from the rather universal value of the constant- 
volume entropy of melting for normal melting elements, 
AS 0.8fc s /atom [H S| (see also |26|, Sec. 22). This 
universality implies the presence of a numerically domi- 
nant and uniform class of potential energy valleys, since 
a distribution of potential energy properties could not 
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be expected to produce the common entropy of melting 
[24j |. The dominant uniform valleys are called random 
valleys. The atomic motion is decomposed into vibra- 
tions in the random valleys, plus complicated but less 
important transit motion, which carries the system be- 
tween valleys. This was argued from the nearly pure 
vibrational values of the atomic motion contribution to 
experimental specific heat and entropy [24l.l25l]. These ar- 
guments are still valid, and still rationalize experimental 
data to an accuracy of a few percent. 



In the original formulation, the transit partition func- 
tion was set to Af r , a constant representing the total num- 
ber of random valleys [24|, [H[ . This expresses the idea 
that all random valleys are equally accessible, and that 
every point in 3V-dimensional configuration space be- 
longs to one and only one potential energy valley [24j . 
However, this approximation for the transit partition 
function is not amenable to systematic improvement in 
order to develop a more accurate theory. With the goal 
of improving transit theory, the available experimental 
high-temperature entropy data for elemental liquids were 
analyzed [27j- The result is a common curve of tran- 
sit entropy for the liquids a naly zed, which represents 
experiment to high accuracy [27|. The purpose of this 
manuscript is to construct a statistical mechanics model 
for the transit partition function and to calibrate the 
model to the experimental transit entropy. The model 
will then provide consistent calibrated equations for the 
transit free energy and all other transit thermodynamic 
functions. 



In Sec. Ill Al the V-T propositions are described. In 
Sec. Ill Bl the V-T Hamiltonian is placed in the framework 
of the formally exact Hamiltonian of a condensed matter 
system ([26j], Sec. 4). In Sec. IIII Al the partition func- 
tion is written, and equations for the internal energy and 
entropy are derived in classical statistical mechanics. In 
Sec. IIII SI the mechanical system representing the liquid 
is identified as the system whose potential energy surface 
consists of the random valleys. The statistical mechanical 
properties of this system which must be accounted for in 
the transit partition function are described. In Sec. IIII Cl 
a model for the transit partition function is constructed, 
and is calibrated to the experimental transit entropy [27| . 
In Sec. IIV Al properties of the calibrated transit model 
are discussed, and in Sec. IIV Bl the model's ability to pre- 
dict experimental internal energy data is shown for liquid 
Na. The physical meaning of the two calibrated model 
parameters is also discussed. Sec. [V] discusses implica- 
tions of two aspects of the present theory. First, with 
the transit model developed here, liquid thermodynamic 
properties can be calculated without adjustable param- 
eters to high accuracy. Second, the symmetry classifica- 
tion of potential energy valleys is expected to be relevant 
beyond the liquid phase, and beyond monatomic systems. 



II. FORMULATION OF V-T THEORY 

A. Basic propositions 

In condensed matter theory, it is generally agreed that 
the potential energy surface is composed of intersecting 
many-atom potential energy valleys. Our hypothesis di- 
vides these valleys into two classes, random and symmet- 
ric, with the following properties in the thermodynamic 
limit (N->oo) dlH. 

(i) The random valleys are macroscopically uniform 
and numerically dominant. Uniformity means for 
any macroscopic dynamical variable, the statistical 
mechanical average is the same for every random 
valley. Numerically dominant means the liquid sta- 
tistical mechanics at T > T m is given entirely by 
the random valleys. 

(ii) The symmetric valleys have a wide range of po- 
tential energy properties. Qualitatively, potential 
properties of symmetric valleys range from crystal 
to liquid values [28l - [3l"j |. Symmetric valleys include 
those with microcrystalline structures; the single 
crystals are also included, and one of them has spe- 
cial status as the ground state structure. 

Current evidence supporting the symmetry classifica- 
tion will be summarized in Sec. El 

Now, since the vibrational motion appears to dominate 
thermodynamic properties of monatomic liquids, we shall 
make this motion the leading term in our Hamiltonian. 
For this purpose, we shall modify the vibrational motion 
to produce a tractable form. Correction for this modifi- 
cation will then become part of the small but complicated 
transit Hamiltonian. 

In a given random valley, at small displacements of the 
atoms from equilibrium, the system potential is quadratic 
in displacements. We define the extended random val- 
ley as the extension to infinity of its harmonic potential 
surface. The motion of atoms in an extended random 
valley is normal-mode vibrational motion, tractable in 
quantum and classical mechanics. Because the extended 
random valleys are uniform in N — ¥ 00, a single such val- 
ley suffices for statistical mechanical calculations. This 
property greatly simplifies the liquid dynamics theory. 

So far we have followed the original formulation. In 
preparation for a general treatment of transits, we shall 
make an accounting of the complete Hamiltonian for a 
condensed matter system ((26[, Sec. 4). 

B. Hamiltonian 

The mechanical system has N atoms in a volume V, 
with periodic boundary conditions on the atomic motion. 
This motion is described by the Hamiltonian "H, where 

u = $ l +n vib +Ht I . (i) 
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3>q(V) is the structural potential of random valleys in 
N — > oo. H v ib describes the vibrational motion in one 
(any) extended random valley, 

^ = £(4 + ^4 (2) 

Here M is the atomic mass, q\ and p\ are respectively 
the normal mode coordinates and momenta, and uj\ are 
the normal mode frequencies, for A = 1, . . . , 3iV. Volume 
dependence of H V ib is contained in the uj\(V). By defi- 
nition of the structure, w\ > for all A, except for the 
three translational modes, for which u>\ = to numerical 
accuracy. 

Utr is the transit Hamiltonian. A transit occurs when 
the system crosses the boundary between two potential 
energy valleys. The motion involves a small local group 
of atoms, and in equilibrium at T > T m , transits are 
occurring at a high rate throughout the liquid. The po- 
tential surface where the system moves in a transit, a 
transit surface, differs locally from the extended random 
valley potential surface. % tr is supposed to express this 
difference for all possible transit surfaces. Our ultimate 
goal is to construct an explicit potential energy function 
for %tv We shall not be able to do that here, but we shall 
be able to construct a simple statistical model that ac- 
counts for what is known about the transit contribution 
to thermodynamics. 

In addition to the Hamiltonian contributions in 
Eq. ([T]), there are three terms which need not be con- 
sidered here. First is H e i, which expresses excitation of 
electrons from their ground state ([26|, Sec. 3). The elec- 
tronic ground state itself is the adiabatic potential for 
atomic motion, and is contained in the three terms on 
the right of Eq. The excitation H e i is determined by 
the electronic density of states evaluated for one (any) 
random structure [32j . This term is important for liq- 
uid metals, contributing 1% — 10% of the liquid internal 
energy and entropy. This term is not explicitly included 
here, because we are discussing only the atomic-motion 
component of liquid dynamics. However, experimental 
information presented here contains proper accounting 
of electronic excitations. 

The next Hamiltonian contribution accounts for the 
interaction between atomic motion and electronic ex- 
citation ([HI, Sec. 4). From calculations for several 
metal crystals [33, 34], we estimate for metallic liquids at 
T > T m that the adiabatic contribution dominates the 
nonadiabatic, and the adiabatic contribution to internal 
energy and entropy is on the order of the experimental 
error in these quantities. We therefore neglect this term, 
because of its smallness. 

Finally, in the actual random valley potential en- 
ergy, there is anharmonicity not associated with transits. 
This is vibrational anharmonicity. In its contribution 
to thermodynamic functions, vibrational anharmonicity 
appears to be a much smaller effect than transits. We 
therefore neglect vibrational anharmonicity, in the sense 



that we shall not attempt to model it explicitly. Ulti- 
mately, however, our calibration of the transit partition 
function to experimental data will include vibrational an- 
harmonicity. The relative effect of this inclusion is pre- 
sumed small (35j . 

A note on the practical calibration of the Hamiltonian 
is useful. The structural and vibrational parts of the 
Hamiltonian are calibrated with the parameters &q(V) 
and {uj\(V)}, respectively. These parameters can be 
calculated from model interatomic potentials [3(1 HH or 
from ab initio electronic structure calculations [3q] . Since 
the calculations are done for finite systems, the results 
have finite- TV errors. It is possible in principle to estimate 
finite- N errors in the Hamiltonian calibration [3(1 HH . 

III. STATISTICAL MECHANICS 

A. Partition function 

Let us write the partition function Z corresponding to 
Eq. (TTJ) for H: 

Z{V,T) = e-KWZ vih (V,T)Z tI (V,T). (3) 
e-WOO i s 

the structural partition function. 
Z v ib(V, T) corresponds to "H v ib, Eq. @, and is fully 
quantum and possesses a classical limit. Z tr (V,T) 
corresponds to 'Htr, and is the primary subject of this 
work. Remaining factors in Z(V,T), not written in 
Eq. ([3]), represent the three remaining Hamiltonian 
contributions mentioned in Sec. Ill B[ 

For simplicity, we shall treat the atomic motion by clas- 
sical statistical mechanics. This is quite accurate for most 
monatomic liquids at T > T^, and the small quantum 
corrections can be estimated ([2q|. Sees. 9 and 17). Since 
the vibrational modes are orthogonal, Z^ is a product 
of single normal-mode functions: 

The factors containing h express kinetic energy. When 
the integrals are done, the result is 
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(5) 



The second equality expresses the definition 1ii(A;b#q) = 
(ln(Hui\)), where ( • • • ) is the average over the set {uj\}. 

Complete expressions for the internal energy U and the 
entropy S are then 

U{V,T) = $ l (V) + 3Nk B T + U tI (V,T), (6) 
S(V,T) = 3Nk B {]n[T/9 l (V)] + l} + 5 tr (V,T).(7) 

The primary theoretical quantities needed to evaluate 
these equations are & (V) for the energy and l (V) for 
the entropy. Additional moments of {lu\} are needed for 
the vibrational quantum corrections. 
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FIG. 1: (Color online) MD data for the mean potential energy 
($) — $o for the random valley system in Na as a function 
of temperature (symbols). The volume is fixed at the vol- 
ume of the liquid at melt at zero pressure. The line is the 
mean vibrational potential energy of (3/2)fc_gT per atom; the 
difference is the energy due to transits. 

B. Mechanical system of the liquid 

Let us define the random valley system as the me- 
chanical system whose potential energy surface consists 
only of random valleys, and all of them. In V-T theory, 
the random valleys make up the configuration space of 
the liquid. The random valley system therefore repre- 
sents the liquid, in the sense that statistical mechanics 
of the random valley system correctly describes the liq- 
uid state. This property holds at all temperatures, at 
T > T m where the liquid is thermodynamically stable, 
and at T < T m where the liquid is metastable with re- 
spect to the crystal. Statistical mechanics of metastable 
states is discussed in (2(|, Sec. 27. 

MD data for temperature dependence of the mean po- 
tential energy, ($) — $q, for the random valley system 
in Na is shown in Fig. [TJ The points are from Fig. 4 of 
[28| . The data are at the fixed volume of the liquid 
at melt at zero pressure, and the notation of volume de- 
pendence will be suppressed. The zero-pressure melting 
temperature of Na is 371.0 K. 

The mean vibrational potential energy is (3/2)fesT per 
atom, and is also graphed in Fig. [TJ The transit internal 
energy per particle, Un, is the difference of the two curves 
in Fig. [JJ 

U tI = «#) - $<) - ^Nk B T. (8) 

Notice this is the potential energy component of Eq. ([6]) . 
We shall not attempt to use this information to calibrate 
a model of Z tI for Na. Rather, we shall use the infor- 
mation to characterize a statistical mechanical model for 
monatomic liquids in general. 

The following significant properties of the random val- 
ley system were observed in [28| . and are reflected in 



Fig.m 

(a) At low temperatures, the equilibrium MD system 
is observed to remain in a single random valley for 
a very long time. The mean potential energy corre- 
sponds to vibrational motion, with anharmonicity 
too small to measure. Equation ([8]) implies [/ tr = 0. 

(b) Upon warming the MD system, a narrow temper- 
ature range is reached where Ut r is observed to in- 
crease from zero. In the same temperature range, 
the self diffusion coefficient D increases from zero. 
Since self diffusion results entirely from transits, 
the coincident appearance of nonzero D and Ut T 
confirms [7 tr as due to transits. 

(c) The high-temperature states in Fig. [TJ are at T > 
T m . The MD values of ($) for these states agree 
with experimental data for liquid Na. 

Let us interpret the preceding observations in terms 
of the statistical mechanics of the random valley system, 
and include the transit entropy in the discussion. This 
will provide the list of transit properties we must account 
for in a statistical mechanics model for Z tT . 

At low temperatures, the random valley system be- 
comes trapped in a single random valley, where the mo- 
tion is entirely vibrational. Since the transit surfaces are 
not visited by the system, both Ut Y and Str are zero. 
Upon warming, the transit rate increases, and at some 
(qualitative) temperature, both Ut T and St r increase from 
zero. This process is thermally activated. The process 
should not be thought of as melting, nor its inverse as 
freezing. Since the random valleys all have the same po- 
tential energy parameters, no first-order phase transition 
is present. With a further increase of temperature, Ut r 
and Str saturate and then begin to decrease. This be- 
havior is not shown in Fig. [TJ but is present in the exper- 
imental internal energy and entropy data of monatomic 
liquids (Sec. E| . 

C. Model for the transit partition function 

The Hamiltonian parameters are U-dependent, while 
T dependence is contained in /3 = (ksT)^ 1 . Volume 
dependence will remain suppressed. Each normal-mode 
configuration integral, Q\, has vibrational contribution 
given by the integral in Eq. (Q}. To include transits, a 
separate transit-surface segment is added to the integral, 
so that 

Qx(T) = J exp (-IpMultfA dq x 

+ |y + ^Jexp(-I3ex)dqx. (9) 

The transit surfaces are b\ < \q\\ < c\, where < b\ < 
c\, and the transit surfaces have potential energy e\ > 0. 
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Evaluation of Eq. © gives 



Qx{T) 



i2nk B T 
Muj\ 




T:k B T 



(10) 

where d\ = c\ — b\. The factor outside the brackets is 
the vibrational contribution, which must be taken with 
the kinetic energy contribution to produce Zvib, Eq. ([5]). 
Then the bracket in Eq. (|10[) is the single normal-mode 
contribution to Z tI , so the total Z tr is then 



Ztr(T) = n 




2Mui\d\ 

exp -/?e A 

TTKb-L 



(11) 



We shall next simplify the model to a form appropriate 
to the work at hand. 

In Eq. pip , each normal mode has two energy param- 
eters, namely e\ and Muj\d\. We shall set each of these 
the same for every mode. It is convenient to use only a 
single energy parameter e, as follows: 



-Mculdl 



(12) 



where /i > is a dimensionless parameter. Equation (jlip 
becomes 



Ztr{T) = [1 

h(T) = M 



rKT)} 
(3e exp( 



-fie). 



(13) 
(14) 



The T dependence of Z tr (T) is entirely contained in /3e. 
This is the source of temperature scaling in the transit 
properties. 

The thermodynamic functions now follow. The 
Helmholtz free energy is _F tr , the internal energy and 
entropy are respectively JTtr and St r , and the constant- 
volume specific heat is Ct r . These are given by 



Ftr = 

U tI = 3Nk B T 



Wk B T1n[l + h(T)] 
(f5e-\)h{T) 



1 + h(T) 



S tT = 3Nk f 



ln[l + h(T)]-\ ^ 2 



\)h{T) 



Ctr = iNk B h(T) 



We 



1)2 



1 + h(T) 

2 



{l + h(T)) 2 l + h(T) 



(15) 
(16) 

(17) 

(18) 



Analysis of the available high-temperature experimen- 
tal entropy data for monatomic liquids has yielded a com- 
mon curve for S tI as a function of Tljh T , where # t r is a 
scaling temperature for each liquid [27| . The common 
curve has a maximum value of Su — 0.8Nk B , located at 
T = 9 tI . Since 9 tr is a material parameter, we must cal- 
ibrate for each liquid independently. Since the T depen- 
dence of Str(T) is contained in j3e, according to Eqs. (TH|) 
and (fTfjl . the model S t r{T) will scale with T/8 tl - if we set 



where v is a dimensionless parameter. Calibration will 
require the determination of /i and v for each liquid. De- 
note by x t ne maximum value of Str(T), which occurs at 
T = #tr, or j3e — v. Then the value and slope of St r (#tr) 
are calibrated to 



Str^tr) = S t i(jl,v) = X, 
Ctr(0tr) = C tI ((l,v)=0. 



(20) 



So fi and v are functions of x only, independent of #tr- If 
X is the same for all liquids, then /i and v are the same 
as well. We set x = 0.8Nk B for monatomic liquids, and 
find 



H = 0.53221, 
v = 1.26452. 



(21) 



Ultimately, the statistical mechanics model is quite 
simple. The complete potential energy surface is modeled 
by two components: vibrational and transit. That the 
vibrational surface lies below the transit surface makes 
Utr and Str go to zero at low temperatures [e > in 
Eqs. (|14|) - (|T7j) ] . That the vibrational surface continues 
above the transit surface causes Ut r and Str to saturate 
and then decrease with increasing temperature [the tran- 
sit term in Eq. pH|) ]. In this way the model possesses the 
required characteristics, as listed in the last paragraph of 

sec. inns 



IV. DISCUSSION OF THE CALIBRATED 
TRANSIT MODEL 

A. Transit entropy 

Figure [2] shows the experimental transit entropy data 
for ten elemental liquids at the fixed volume (the 
volume of the liquid at melt at zero pressure), as a func- 
tion of T/8tr- Also shown is St T (T/9 t r) from the sta- 
tistical mechanics model, Eq. (|17[) . with the calibration 
Eq. (|2ip . The agreement of model with experiment is 
excellent; scatter of the points from the line is of order 
0.01 fcs/atom, or around 0.1% of the total entropy. How- 
ever, we do not suppose the model is entirely "correct" 
to this level of accuracy. Systematic errors of the experi- 
mental entropy data are larger than the scatter in Fig. [5J 
and these errors can affect the shape of the St T (T) curve. 
This uncertainty is reflected in the estimate of up to 10% 
error in the fitted values of 6*tr (27| . 

To understand Fig. [5J one needs to know how it is made 
[27j . The zero-pressure experimental data S cxp t(V,T) 
are corrected, by means of additional experimental data, 
to the fixed volume V^. This produces a data set for 
Sexpt{VjniT) for each liquid. This data set is then fitted 
to the right side of Eq. ([7]), with l a variable parameter. 
In the fitting process, the maximum of each S , tr (U T ^,T) 
curve is fixed at 0.8A:B/atom, and the temperature at 
the maximum is denoted 9 tr . The fitted characteristic 



vk 



BVtr, 



(19) temperatures 9 and 9 tl - are listed in [27 
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FIG. 2: (Color online) Experimental transit entropy data for 
ten elemental liquids as a function of temperature scaled by 
the characteristic transit temperature (9 tr . The data are from 
Fig. 1 of [27]]. The solid line is the transit model presented 
here, Eq. (fT7|l . calibrated as indicated in Eq. ([21]). 
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FIG. 3: (Color online) Transit energy per atom for Na, deter- 
mined from MD calculations and experiment, as a function 
of temperature at fixed volume. The solid line is the tran- 
sit model presented here, Eq. (|16[l . calibrated as indicated in 
Eq. (21). 



The presence of temperature scaling in the statistical 
mechanics model is informative. It results from the sim- 
plicity of the model, in that T appears in only one form in 
Z tT , namely in the form j3e [see Eqs. ([13]) and (|TI|)]. The 
previously-discovered scaling temperature #tr is then in- 
troduced via Eq. (|19p , leaving v as the parameter in place 
of e. The following dimensionless functions are then func- 
tions only of T/6t t : @e, h, Str/fes, Ctr/fce, F tr /fcsT, and 
Utr/ksT. This scaling behavior should prove useful in 
studying the transit thermodynamic functions, especially 
since they are relatively small contributions. 

Finally in Fig. [5] we note that the model rises slightly 
above the data at T > 26 tI . This discrepancy can be due 
to experimental errors in the data analysis. However, 
the discrepancy might be significant, and due to transits 
at very high potential energy. Here, the major transit 
correction is the removal of the vibrational surface that 
was extended beyond the random valley boundary. That 
correction is modeled in a study of liquid Hg at high tem- 
peratures [37j ■ This effect would break the T/9 tT scaling 
of Str- But the effect seen in Fig.[5]is too small to justify 
modeling. 

B. Transit internal energy 

Experimental data for the zero-pressure internal en- 
ergy U (V, T) for liquid Na are corrected to the fixed vol- 
ume Vjn, and the experimental [/ t r(^, T) is extracted by 
means of Eq. ©. For this, the value $ (V^) = 0.33±0.05 
mRy/atom is estimated from early MD data dHHl]. & 
is measured relative to the thermodynamic zero of energy, 
which is the energy of the crystal at zero temperature 
and pressure. The experimental data for UniV^jT) are 
graphed in Fig. [3] MD data for U tl (V^,T) are obtained 
from Fig. Q] by means of Eq. (|5J), and are also graphed 



in Fig. [3] The model curve is calculated from Eq s. ([16]) 
and ([21]), together with 9 tI = 570 K for liquid Na [27|. 

In Fig. [3] the model curve rises above experiment at 
T > 800 K. This is due to error in the volume correction 
of the experimental energy. That volume correction is 
extremely difficult to evaluate, hence becomes inaccurate 
at a relatively low temperature. Also in Fig. [3] at T m , 
the MD curve is lower than experiment by 1.4% of the 
total experimental internal energy. This is likely due to 
error in the Na interatomic potential, plus a few smaller 
contributing errors. Overall, the discrepancies in Fig. [3] 
are remarkably small, and can be ignored in our analysis. 

Let us first compare the model with experiment in 
Fig. [3] Since the model agrees with experiment for S tI (T) 
at T > T m , Fig. [2] it must also agree with experiment 
for Utr(T) at T > T m , with the possible exception of an 
error in a constant of integration. The constant of inte- 
gration can be taken as the value of Ut T (T m ). This being 
given correctly by the model constitutes a verification of 
the model, independent of the verification provided by 
the entropy in Fig. [2] 

Let us next compare the model with MD in Fig. [3] 
The comparison is proper since both are classical. One 
immediately sees that /7t r (T) increases from zero at the 
same temperature for both the model and MD. Notice 
also that S*t r (T) will increase from zero at the same tem- 
perature as Utr(T) does, both for MD data and for the 
model. In Fig. [3] agreement of the model with MD at 
T < T m is not entirely independent of the agreement with 
experiment at T > T m . It is nevertheless quite satisfac- 
tory that the transit model, calibrated at T > T m , goes 
to zero at the correct temperature (for Na) well below 
T 

In Sec. ITHBl it was argued that S tr (T) and U tr (T) 
must vanish at temperatures where transits are not ther- 
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mally activated in the random valley system. Hence, 
S tr (T) contains no additive constant. This implies that 
the random valley multiplicity does not contribute to 
Str(T), so that the effective number Af r of random valleys 
satisfies \nj\f r < O(N). The agreement between theory 
and experiment in Figs. [2] and [3] provides support for this 
conclusion. 

The model parameters have physical meaning, which 
is sharpened by the calibration process. Our tentative 
presumption is that this discussion applies to monatomic 
liquids in general. The potential energy of each effec- 
tive transit surface is e, measured from the structural 
potential $o, the same level from which vibrational en- 
ergy is measured. From the calibration of v, Eq. (f2"Tj) . we 
have e = 1.265fc,e#tr. Hence the transit surface is easily 
accessible at liquid temperatures, /i is related to the ef- 
fective length of the transit surface in each q\ direction, 
and [i fixes the magnitude of h(T), Eq. (|14l) . The ra- 
tio of transit to vibrational contributions in the internal 
energy and entropy is of order h{T), which is small com- 
pared to 1 at all T. The ratio can be calculated exactly 
from the equations of Sec. IIIII This provides a quanti- 
tative statement of our introductory message, that the 
transit contribution to liquid thermodynamic properties 
is small, but important for accurate work. 



V. THEORETICAL PREDICTIONS AND 
VERIFICATIONS 

With the calibrated transit model of Sec. IIIII it is pos- 
sible to evaluate thermodynamic properties of elemen- 
tal liquids without adjustable parameters. We consider 
specifically the internal energy and entropy, whose V-T 
theory formulas are given in Eqs. (JSJ) and ([7]). To date, 
three independent tests verify these formulas to an ac- 
curacy within experimental error, and are therefore con- 
sistent with zero theoretical error. By comparing ab ini- 
tio calculations of 9 l with the values determined by fit- 
ting the experimental entropy, it is shown that Eq. ([7]) 
agrees with experiment for Na and Cu (27j . Equation 
©, with MD evaluation of and with U tI (T) from 
Eq. (|I6p. agrees with experiment for the internal energy 
of Na (Fig. 13]). Continued testing in this way will reveal 
the overall accuracy of the theory, and will uncover any 
cases where the theory needs significant correction. 

In comparing theory and experiment for the liquid en- 
tropy, the temperature dependence of experimental en- 
tropy is already accounted for by the transit model for 
the ten liquids in Fig. [2] For other liquids, tempera- 
ture dependence of experimental entropy is independent 
information, so that the entire temperature dependence 
of the internal energy or entropy will test the theory. 
In the analysis of experimental data reported in this 



paper, electronic excitation contributions are evaluated 
from free electron theory. A more accurate calculation 
is based on the electronic density of states evaluated for 
a random structure (Sec. Ill Bp . The random structure 
density of states will be necessary for transition metals, 
where electronic excitation contributions are much larger 
than in the nearly-free-electron metals. This provides 
an additional theoretical prediction which can be tested. 
The same theory was used to isolate the anharmonic vi- 
brational contribution to entropy in the transition metal 
crystals (40| . 

Our development of monatomic liquid dynamics the- 
ory is strongly based on the symmetry classification of 
potential energy valleys, Sec. Ill Al The same random 
valleys important for liquid theory will continue to domi- 
nate at higher temperatures, where the system undergoes 
the broad liquid-to-gas transition [4l|. Moreover, all the 
potential valleys, random and symmetric, will contribute 
to the description of amorphous solids at temperatures 
well below T m . In view of this extended application of 
the symmetry classification hypothesis, it is worthwhile 
to investigate to what classes of materials it applies. 

The symmetry classification hypothesis is now well ver- 
ified for Na at the density of the liquid at melt [13, [H| ■ 
The complete distribution of structural potentials $0 
shows the narrow but dominant random peak, the broad 
symmetric distribution, and the crystal 31]. Other mea- 
sures can distinguish randoms from symmetries, e.g., 
Voronoi analysis and pair correlations [29) . Work in 
progress shows the dominance and uniformity of ran- 
dom structures for Al and Cu. The logical next step 
is to test the hypothesis for monatomic systems in gen- 
eral. The characteristic random structures are expected 
to vary with volume for one element, and to vary from 
one element to another. 

In principle, the symmetry classification of potential 
energy valleys should also apply to more complicated sys- 
tems, e.g., to alloys, compounds, and molecular systems. 
Evidence has been cited for the presence of a dominant 
and uniform class of potential valleys, i.e., the random 
class, in a variety of MD systems (see (30| and references 
[56-62] quoted there) . If the random valley class is found 
to be present in complex liquids, we shall have moved 
a step closer to a Hamiltonian formulation for such sys- 
tems. 
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